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Abstract. We study the two primary families of spaces of finite element differential forms with respect 
to a simplicial mesh in any number of space dimensions. These spaces are generalizations of the classical 
finite clement spaces for vector fields, frequently referred to as Raviart— Thomas, Brczzi— Douglas-Marini, and 
Nedelec spaces. In the present paper, we derive geometric decompositions of these spaces which lead directly 
to explicit local bases for them, generalizing the Bernstein basis for ordinary Lagrange finite elements. The 
approach applies to both families of finite element spaces, for arbitrary polynomial degree, arbitrary order 
of the differential forms, and an arbitrary simplicial triangulation in any number of space dimensions. A 
prominent role in the construction is played by the notion of a consistent family of extension operators, 
which expresses in an abstract framework a sufficient condition for deriving a geometric decomposition of a 
finite element space leading to a local basis. 



1. Introduction 

The study of finite element exterior calculus has given increased insight into the construction of stable 
and accurate finite clement methods for problems appearing in various applications, ranging from electro- 
magnetics to elasticity. Instead of considering the design of discrete methods for each particular problem 
separately, it has proved beneficial to simultaneously study approximations of a family of problems, tied 
together by a common differential complex. 

To be more specific, let C R" and let HA k (Q) be the space of differential k forms uo on ft, which is in 
L 2 , and where its exterior derivative, duj, is also in L 2 . The L 2 version of the de Rham complex then takes 
the form 

-» HA°(fl) i HA 1 (ft) ■ ■ ■ -i HA n (fl) -» 0. 
The basic construction in finite element exterior calculus is of a corresponding subcomplex 
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where the spaces A\ are finite dimensional subspaces of HA k (Q) consisting of piecewise polynomial differen- 
tial forms with respect to a partition of the domain £1. In the theoretical analysis of the stability of numerical 
methods constructed from this discrete complex, bounded projections Hh ■ HA k (fl) — > Ajj are utilized, such 
that the diagram 
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commutes. For a general reference to finite clement exterior calculus, we refer to the survey paper [2J, and 
references given therein. As is shown there, the spaces are taken from two main families. Either A^ is one 
of the spaces V r A k (T) consisting of all elements of HA k (il) which restrict to polynomial k- forms of degree 
at most r on each simplex T in the partition T, or A^ — r P~A k (T), which is a space which sits between 
P r A k (T) and V r -iA k (T) (the exact definition will be recalled below). These spaces are generalizations of 
the Raviart-Thomas and Brezzi-Douglas-Marini spaces used to discretize -ff(div) and if(rot) in two space 
dimensions and the Nedelec edge and face spaces of the first and second kind, used to discretize ff (curl) and 
i/(div) in three space dimensions. 

A key aim of the present paper is to explicitly construct geometric decompositions of the spaces V r A k (T) 
and V~A k (T) for arbitrary values of r > 1 and k > 0, and an arbitrary simplicial partition T of a polyhedral 
domain in an arbitrary number of space dimensions. More precisely, we will decompose the space into 
a direct sum with summands indexed by the faces of the mesh (of arbitrary dimension), such that the 
summand associated to a face is the image under an explicit extension operator of a finite-dimensional space 
of differential forms on the face. Such a decomposition is necessary for an efficient implementation of the 
finite element method, since it allows an assembly process that leads to local bases for the finite element 
space. The construction of explicit local bases is the other key aim of this work. 

The construction given here leads to a generalization of the so-called Bernstein basis for ordinary poly- 
nomials, i.e., 0-forms on a simplex T in K™, and the corresponding finite element spaces, the Lagrange finite 
elements. See Section pO] below. This polynomial basis is a well known and useful theoretical tool both in 
finite element analysis and computational geometry. For low order piecewise polynomial spaces, it can be 
used directly as a computational basis, while for polynomials of higher order, this basis can be used as a 
starting point to construct a basis with improved conditioning or other desired properties. The same will be 
true for the corresponding bases for spaces of piecewise polynomial differential forms studied in this paper. 

This paper continues the development of geometric decompositions begun in [2j Section 4]. In the present 
paper, we give a prominent place to the notion of a consistent family of extension operators, and show that 
such a family leads to a direct sum decomposition of the piecewise polynomial space of differential forms 
with proper interelement continuity. The explicit notion of a consistent family of extension operators is new 
to this paper. We also take a more geometric and coordinate-independent approach in this paper than in 
[2], and so are able to give a purely geometric characterization of the decompositions obtained here. The 
geometric decomposition we present for the spaces V~ A k here turns out to be the same as obtained in [2J, 
but the decomposition of the spaces V r A k obtained here is new. It improves upon the one obtained in [2J, 
since it no longer depends on a particular choice of ordering of the vertices of the simplex T, and leads to a 
more canonical for V r A k . 

The construction of implementable bases for some of the spaces we consider here has been considered 
previously by a number of authors. Closest to the present paper is the work of Gopalakrishnan, Garcfa- 
Castillo, and Demkowicz [4]. They give a basis in barycentric coordinates for the space V~A , where T is 
a simplex in any number of space dimensions. In this particular case, their basis is the same as we present 
in S ection [9} In fact, Table 3.1 of [I] is the same, up to a change in notation, as the left portion of Table 



9.2| of this paper. As will be seen below, explicit bases for the complete polynomial spaces V r A k are more 
complicated than for the V~ A k spaces. To our knowledge, the basis we present here for the V r A k spaces 
have not previously appeared in the literature, even in two dimensions or for small values of r. 

Other authors have focused on the construction of p-hierarchical bases for some of the spaces considered 
here. We particularly note the work of Ainsworth-Coyle [Q, Hiptmair [5], and Webb [7J. In [T], the authors 
construct hierarchical bases of arbitrary polynomial order for the spaces we denote V r A k , k = 0, . . . , 3, r > 1, 
and T a simplex in three dimensions. In section 5 of [5], Hiptmair considers hierarchical bases of V~ A k for 
general r, k, and simplex dimension. In 7, Webb constructs hierarchical bases for both V r A k and V~A k , 
for k = 0, 1 in one, two, and three space dimensions. The approaches of these three sets of authors differ. 
Even when adapted to the simple case of zero-forms, i.e., Lagrange finite elements, they produce different 
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hierarchical bases, from among the many that have been proposed. Our approach is quite distinct from these 
in that we are not trying to find hierarchical bases, but rather we generalize the explicit Bernstein basis to 
the full range of spaces V,A k and V~A k . 

In the present work, by treating the V r A k and V~ A k families together, and adopting the framework of 
differential forms, we are able to give a presentation that shows the close connection of these two families, 
and is valid for all order polynomials and all order differential forms in arbitrary space dimensions. Moreover, 
the viewpoint of this paper is that the construction of basis functions is a straightforward consequence of 
the geometric decomposition of the finite clement spaces, which is the key ingredient needed to construct 
spaces with the proper intcr-clcment continuity. Thus, the main results of the paper focus on these geometric 
decompositions . 

An outline of the paper is as follows. In the next section, we define our notation and review material we will 
need about barycentric coordinates, the Bernstein basis, differential forms, and simplicial triangulations. The 
V r A k and V~ A k families of polynomial and piecewise polynomial differential forms are described in Section[3] 
In Section [4] we introduce the concept of a consistent family of extension operators and use it to construct 
a geometric decomposition of a finite element space in an abstract setting. In addition to the Bernstein 
decomposition, a second familiar decomposition which fits this framework is the dual decomposition, briefly 
discussed in Section [5] Barycentric spanning sets and bases for the spaces V r A k (T) and V~A k jT) and the 
corresponding subspaces V r A k (T) and T~ A k {T) with vanishing trace are presented in Sectional The main 
results of this paper, the geometric decompositions and local bases, are derived in Sections [7] and [8] for 
T>-A k (T) and V r A k (T). Finally, in Section [9] we discuss how these results can be used to obtain explicit 
local bases, and tabulate such bases in the cases of 2 and 3 space dimensions and polynomial degree at most 
3. 



2. Notation and Preliminaries 



2.1. Increasing sequences and multi-indices. We will frequently use increasing sequences, or increasing 
maps from integers to integers, to index differential forms. For integers j, k, I, m, with < k — j < m — /, we 
will use S(j : k,l : m) to denote the set of increasing maps {j, . . . , k} — ► {I, . . . , m}, i.e., 

E(j : k, I : m) = { o : {j, . . . , k} -> {I, . . . , m} | a(j) < a(j + 1) < • • • < a(k) }. 

Furthermore, [cr] will denote the range of such maps, i.e., for cr G S(j : k,l : m), Jer] = { a{i) | i = j, . . . , k}. 
Most frequently, we will use the sets E(0 : k, : n) and E(l : k, : n) with cardinality or ), 

respectively. Furthermore, if a G E(0 : k, : n), we denote by cr* £ S(l : n — k, : n) the complementary 
map characterized by 

(2.1) HU[(7*]={0,l,...,n}. 

On the other hand, if cr G S(l : k, : n), then cr* G S(0 : n — k, : n) is the complementary map such that 
fl2Tj ) holds. 

We will use the multi-index notation a G Nq , meaning a = [a%, • • • , a n ) with integer a» > 0. We define 
x a = x" 1 ■ ■ - a;"", and \a\ :— J2 a i- We will also use the set Nq" of multi-indices a = (a n , ■ ■ ■ ,a„), with 
n . The support [a] of a multi-index a is { i | a, > }. It is also useful to let 

[a,aj = H U [cr], a G Nq " , a G E(j : k, I : m). 

If Cl C R™ and r > 0, then V r (Q) denotes the set of real valued polynomials defined on of degree 
less than or equal to r. For simplicity, we let Vr = V r (M. n ). Hence, if has nonempty interior, then 
dim'P. r (Sl) = dimVr = ( r ^™)- The case where f2 consists of a single point is allowed: then V r (Q) = K for all 
r > 0. For any ft, when r < 0, we take V r (£l) = {0}. 
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2.2. Simplices and barycentric coordinates. Let T € M. n be an n-simplex with vertices xq, x\, . . . , x n 
in general position. We let A(T) denote all the subsimplices, or faces, of T, while A^T) denotes the set 
of subsimplices of dimension k. Hence, the cardinality of Afc(T) is (^+1)- We will use elements of the set 
E(j : fc, : n) to index the subsimplices of T. For each a € E(j : k, : n), we let /o- € A(T) be the closed 
convex hull of the vertices x a ^, . . . , 2Wfc), which we henceforth denote by [a: CT (j), ■ ■ ■ , £<r(fc)]- Note that there 
is a one-to-one correspondence between Afc(T) and S(0 : fe, : n). In fact, the face f a is uniquely determined 
by the range of a, Jcr]. If / = f a for a e E(j : fc,0 : n), we let the index set associated to / be denoted by 
1(f), i.e., 1(f) = {a}. If / € A fe (T), then /* e A„_ fe _i(T) will denote the subsimplex of T opposite /, i.e., 
the subsimplex whose index set is the complement of 1(f) in { 0, 1, . . . , n }. Note that if a <E E(0 : k, : n) 
and / = f a , then /* = f a ,. 

We denote by Ajf, Xj , . . . , the barycentric coordinate functions with respect to T, so Xf € V\(T) is 
determined by the equations Xf(xj) = Sij, < i,j < n. The functions Xf form a basis for Vi(T), are 
non-negative on T, and sum to 1 identically on T. Moreover, the subsimplices of T correspond to the zero 
sets of the barycentric coordinates, i.e., if / = f a for a E S(0 : k, : n), then / is characterized by 

f = {xeT\Xf(x)=Q, ie{a*}}. 

For a subsimplex / G A(T), the barycentric coordinates functions with respect to /, {X{} ieI (f) C ~Pi(f), 
satisfy 

(2.2) x{ = tr tj Xj , i £ 1(f)- 



Here the trace map tr^,/ : Vi(T) — > V\(f) is the restriction of the function to /. Due to the relation (2.2 1, 
we will sometimes omit the superscript T or /, and simply write Xi instead of Xf or X{ . Note that, by 
linearity, the map X{ — > Af , i € defines a barycentric extension operator Ej T ; V\(f) — * Vi(T), which 

is a right inverse of try,/- The barycentric extension Ej T p can be characterized as the unique extension of 
the linear polynomial p on / to a linear polynomial on T which vanishes on /*. 



2.3. The Bernstein decomposition. Let T = [xq, X\, . . . , x n ] C be as above and {Ai}™ =0 C V\(T) the 
corresponding barycentric coordinates. For r > 1, the Bernstein basis for the space V r (T) consists of all 
monomials of degree r in the variables A^, i.e., the basis functions are given by 

(2.3) { A" = A^A" 1 • • • A"" I a <E Nq", |a|=r}. 

(It is common to take the scaled barycentric monomials (n\/al)X a as the Bernstein basis elements, as in [6], 
but the scaling is not relevant here, and so we use the unsealed monomials.) Of course, for / £ A(T), the 
space V r (f) has the corresponding basis 

{(A'HaGNr, H = r, [a] C 1(f) }. 

Hence, from this Bernstein basis, we also obtain a barycentric extension operator, E = E^ T : V r (f) — > V r (T), 
by simply replacing X{ by in the bases and using linearity. 

We let V r (T) denote the subspace of V r (T) consisting of polynomials which vanish on the boundary of T 
or, equivalently, which arc divisible by the corresponding bubble function Ao • • • A„ on T. Alternatively, we 
have 

(2.4) V r (T) = span{A« | a € N° : ", \a\ = r, [a] = {0, . . . , n} }. 

Note that multiplication by the bubble function establishes an isomorphism V r -n-i(T) — V r (T). 



The Bernstein basis (2.3 1 leads to an explicit geometric decomposition of the space V r (T). Namely, we 
associate to the face /, the subspace of V r (T) that is spanned by the basis functions A Q with [a] = 1(f). 
We then note that this subspace is precisely E[P r (f)], i.e., 

(2.5) E[V f (f)] = span{ A Q | a S \a\ - r, {a} = 1(f) }. 
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Clearly, 

(2.6) V r (T)= E[V r (f)], 

/£ A(T) 

which we refer to as the Bernstein decomposition of the space V r (T). This is an example of a geometric 
decomposition, as discussed in the introduction. An illustration of the decomposition (2.6 1 is given in 
Figure |2Jj 




Figure 2.1. The Bernstein basis of V± (T) for a triangle T. One basis function is associated 
with each vertex, three with each edge, and three with the triangle. The basis functions 
associated with any face / are obtained by extending basis functions for V^f) to the triangle. 



Moreover, the extension operator E may also be characterized geometrically, without recourse to barycen- 
tric coordinates. To obtain such a characterization, we first recall that a smooth function u : T — > K is said 
to vanish to order r at a point x if 

(d a u)(x) = 0, a€N£, |a|<r-l. 

We also say that u vanishes to order r on a set g if it vanishes to order r at each point of g. Note that the 
extension operator E — Ej T has the property that for any /i S V r (f), Efi vanishes to order r on /*. In 
fact, if we set 

V r (T, f) = {ue V r (T)\uj vanishes to order r on /* }, 

we can prove 

Lemma 2.1. V r (T, f) — E[P r (f)] and for fi G Pr(/)j E/i = Ef T fi can be characterized as the unique 
extension of fj, to V r (T,f). 

Proof. It is easy to see that E[P r (f)] C V r (T,f). To establish the reverse inclusion, we observe that if 
/* = {xi}, then UJ S V r {T) vanishes at /* if and only if it can be written in the form 

\cx\—r 
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where the sum is restricted to multi-indices i for which on = 0. For a more general set /*, this fact will be 
true for any i £ and hence 

uj= X] CqAQj 

\a\—r 
I«]CX(/) 

and so w G E[P r (f)]. □ 



We also note that it follows immediately from Lemma 2.1 that the space E[P r (f)] appearing in the 
Bernstein decomposition ( |2.6[ ) is characterized by 

E[Pr(f)] = {u G V r {T) | to vanishes to order r on /*, tr TJ w G V r (f)}. 



In this paper, we will establish results analogous to those of this section for spaces of polynomial differential 
forms, and in particular, direct sum decompositions of these spaces analogous to the Bernstein decomposition 



2.4. Differential forms. Next we indicate the notations we will be using for basic concepts related to 
differential forms. See §2] or the references indicated there for a more detailed treatment. For k > 0, we 
denote by Alt fe V the set of real- valued, alternating k- linear maps on a vector space V (with Alt V = R). 
Hence, Alt fc V is a vector space of dimension ( d ™ y )- The exterior product, or the wedge product, maps 
Alt 3 V x Alt' 1 ' V into Alt 3+k V. If w G Alt fe V and v £ V, then the contraction of w with v, lujv £ Alt fe_1 V, 
is given by ujv(v\,. . . ,v k -i) = u(v,vx, . . .,v k -i). 

If is a smooth manifold (e.g., an open subset of Euclidean space), a differential fc-form on is a map 
which assigns to each iGSlaii element of Alt fc T x fl, where T x fl is the tangent space to f2 at x. In case / is 
an open subset of an affine subspace of Euclidean, all the tangents spaces T x f may be canonically identified, 
and we simply write them as Tf. 

We denote by A k (£l) the space of all smooth differential fc-forms on fL The exterior derivative d maps 
A fc (f2) to A fc+1 (f2). It satisfies do d = 0, so defines a complex 

-» A {SI) -i A\Q) ±> ■ ■ ■ -i A n (fi) -» 0, 

the de Rham complex. 

If F : SI — > is a smooth map between smooth manifolds, then the pullback F* : A fc (J7') — > A fe (J7) is 
given by 

(F*w)»(«i,«2, . . . ,« fc ) = UFwiDFsfa^DFvfo), . . . , DF x (v k )), 
where the linear map DF X : T x St — » Tp^jfJ' is the derivative of F at x. The pullback commutes with the 
exterior derivative, i.e., 

F*{dLo) = d{F*u), uj£A k (Sl'), 
and distributes with respect to the wedge product: 

f'(wAtj) = F*uj A F*rj. 
We also recall the integral of a fc-form over an orientable /c-dimcnsional manifold is defined, and 

(2.7) [ F*lu= [ u, Lj£A n (Sl'), 

Jn Jw 

when F is an orientation-preserving diffeomorphism. 

If SI' is a submanifold of St, then the pullback of the inclusion SI' SI is the irace map tr^.w '■ A k (St) — > 
A fc (J7'). If the domain 17 is clear from the context, we may write tr^' instead of ti'o.si', and if SI' is the 
boundary of SI, dSl, we just write tr. Note that if SI' is a submanifold of positive codimension and k > 0, 
then the vanishing of Trn.n' w on ^' f° r w G A k (Sl) does noi imply that to x £ Alt k T x Sl vanishes for x £ St', 
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only that it vanishes when applied to fc-tuples of vectors tangent to f2', or, in other words, that the tangential 
part of lu x with respect to T x il' vanishes. 

If O is a subset of R n (or, more generally, a Riemannian manifold), we can define the Hilbert space 
L 2 A fc (f2) z> A fc (f2) of L 2 differential fc-forms, and the Sobolev space 

HA k (n) - {cue L 2 A k (Q)\du e L 2 A k+1 (fl)}. 

The L 2 de Rham complex is the sequence of mappings and spaces given by 

(2.8) HA°(Q) ± HA\n) -i • • • -4 HA n (Q) 0. 

We remark that for f2 C E™, ff A°(f2) is equal to the ordinary Sobolev space i? 1 (fi) and, via the identification 
of Alt" R n with R, HA n (Q) can be identified with L 2 (fi). Furthermore, in the case n = 3, the spaces Alt 1 M 3 
and Alt 2 R 3 can be identified with R 3 , and the complex ( |2.8[ ) may be identified with the complex 

_ ^ (curl; n) ^> H(div; fi) ^> T 2 (ft) -» 0. 

2.5. Simplicial triangulations. Let f2 be a bounded polyhedral domain in R n and T a finite set of n- 
simplices. We will refer to T as a simplicial triangulation of CI if the union of all the elements of T is the 
closure of f2, and the intersection of two is either empty or a common subsimplex of each. For < j < n, 
we let 

n 

A 3 (T) = |J A,(T) and A(T) = |J A,(T). 

TeT j=0 

In the finite element exterior calculus, we employ spaces of differential forms u> which are piecewise smooth 
(usually polynomials) with respect to T, i.e., the restriction co\t is smooth for each TeT. Then for 
/ € Aj(T) with j > k, ttf lu may be multi- valued, in that we can assign a value for each T G T containing 
/ by first restricting u> to T and then taking the trace on /. If all such traces coincide, we say that tr/cj is 
single- valued. The following lemma, a simple consequence of Stokes' theorem, cf. [5] Lemma 5.1], is a key 
result. 

Lemma 2.2. Let ui E L 2 A k (il) be piecewise smooth with respect to the triangulation T. The following 
statements are equivalent: 

(1) to e HA k (fl), 

(2) tifui is single-valued for all f G A„_i(T), 

(3) tr f lu is single-valued for all f G A,-(T), k < j < n— 1. 

As a consequence of this theorem, in order to construct subspaces of HA k (£l), consisting of differential 
forms lu which arc piecewise smooth with respect to the triangulation T, we need to build into the construction 
that tr/ lu is single- valued for each / G Aj (T) for k < j < n — 1 . 



3. Polynomial and piecewise polynomial differential forms 

In this section we formally define the two families of spaces of polynomial differential forms V r A k and 
V~ A k . These polynomial spaces will then be used to define piecewise polynomial differential forms with 
respect to a simplicial triangulation of a bounded polyhedral domain in R n . In fact, as explained in 2, §3.4], 
the two families presented here are nearly the only affine invariant spaces of polynomial differential forms. 
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3.1. The space V r A k '. Let ft be a subset of W\ For < k < n, we let V r A k (n) be the subspace of A fe (ft) 
consisting of all uj G A fe (ft) such that io{v\, V2, ■ ■ ■ , v/~) G 7-V(ft) for each choice of vectors V\,v%, . . .v^ G R™. 
Frequently, we will write V r A k instead of V r A k (M. n ). The space V r A k is isomorphic to V r ® Alt fe and 

'r + n\ / n\ / r + k\ ( n + r x 
n J \k J \ r J \n — k j 

Furthermore, if ft C K™ with nonempty interior, then dim7 : VA fc (ft) = dim'P r A fc . 



(3.1) dim V r A fe = dim V r x dim Alt fc TD> " 



If T is a simplex, we define 

A,A fc (T) = { uj e V r A k (T) I tr u = }. 
In the case k = 0, this space simply consists of all the polynomials divisible by the bubble function Ao • • • A„, 
so 

(3.2) V r A°(T)=Vr-n~iA°(T). 
For k = n, the trace map vanishes, so we have 

(3.3) V r A n (T)=V r A n (T). 

3.2. The space V~A k . The Koszul differential k of a differential fc-form 10 on K n is the (fc — l)-form given 
by 

(kuj) x (vx, . . . ,v k -i) = ^(^(a;),^!,.. .,v k -i), 

where X(x) is the vector from the origin to x. For each r, k maps V r -iA k to V r A k ~ 1 , and the Koszul 
complex 

-> V r -nA n ^— > Pr-n+lA"- 1 ^— > • • • ^— ► ^A K 0, 

is exact. Furthermore, the Koszul operator satisfies the Leibniz relation 

(3.4) k(uj A 77) = (/tu) Ar/+ (-l)'wA (k»j), w G A fe , ry e A'. 
We define 

■p.7A fc = P r rA i: (M") = V r -iA k + nV r -iA k+1 . 
From this definition, we easily see that V~A° = ^A and V~A n = V r -\A n . However, if < k < n, then 

V r -iA k C 7 , J TA fe C V r A k . 
An important property of the spaces V~A k is the closure relation 

(3.5) V~A k A V~A l C p- +s A fc+i . 

A key identity relating the Koszul operator k with the exterior derivative d is the homotopy relation 

(3.6) (dn + Kd)u = (r + k)u, uj G H r A k , 
where 7i r A k is the space of homogeneous polynomial fc-forms of degree r. 

Using the homotopy relation and the exactness of the Koszul complex, we can inductively compute the 
dimension of V~A k as 



(3.7) dim?VA fc 



r+fc — l\ / n + r 
k ) \n — k 



If ft C R™, then V r ~A k (Cl) denotes the restriction of functions in V~A k to ft, which implies that the space 
V~ A fe (ft) is isomorphic to V~A k if ft has nonempty interior. Finally, we remark that although the Koszul 
operator k depends on the choice of origin used to associate a point in R™ with a vector, the space V~A k 
is unaffected by the choice of origin. We refer to [5] for more details on the spaces V~A k . In particular, if 
T C R n is a simplex and / G Aj(T) then tx f V-A k {T) = V-A k (f), where the space V-A k (f) = V-A k (W) 
depends on /, but is independent of T. 

For a simplex T, we define 

V~A k (T) = {uje 7VA fe (T) | tru = 0}. 
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From the Hodge star isomorphism, we have that V r -n-iA°{T) = V r - n -iA n (T) = V~_ n A n (T) and that 



V r A n (T) = V r A°(T) = V~A°(T). Therefore (Q and Q become 

(3.8) V r A°(T) = V-_ n A n (T), V r A n (T) Si P"A (T). 
These are the two extreme cases of the relation 

(3.9) r r A k (T)=V^_ n+k A n - k (T), 0<k<n. 
But (|3.8h can also be written 



V-A°(T) = 7V»-iA"(T), V-A n {T) - V r - X A\T), 
(where we have substituted r — 1 for r in the second relation) , which are the extreme cases of 
(3.10) V-A k (T) £ 7V n+fe _iA"- fc (r), < k < n. 



That the isomorphisms in (3.9 1 and (3.101 do indeed exist for all k follows from Corollary 5.2 below 



3.3. The spaces 7-VA fe (T) and V r A k (T). For T a simplicial triangulation of a domain S K™, we define 

V r A k (T) = {w£ L 2 A fc (ft) | w| T e V r A k (T) VT e T, 

tr/ a; is single- valued for / G Aj(T), k < j < n — 1}, 

and define V^A k (T) similarly. In view of Lemma 



2.2 



we have 



V r A k (T) = {uje HA k (tt) | w|t € P r A fe (T) VT e T}, 
V-A k (T) = {u e HA k (n) \ lo\t S V-A k (T) VT e T }. 

4. Consistent extension operators and geometric decompositions 

Let T he & simplicial triangulation of Q, C R™, and let there be given a finite dimensional subspace X(T) 
of A fe (T) for each T e T. In this section we shall define the notion of a consistent family of extension 
operators, and show that it leads to the construction of a geometric decomposition and a local basis of the 
finite element space 

(4.1) X{T) = {cue L 2 A k (n) | oj\ T £ X(T) VT e T, tr f uj is single-valued for / € A(T) }. 

X(T) c HA k (fl). 



We note that as a result of Lemma 



2.2 



For the Lagrange finite element space V r (T) = 7 , r A°(T), both the Bernstein basis discussed in Section 2.3 
and the dual basis discussed in the next section arise from this construction. One of the main goals of this 
paper is to generalize these bases to the two families of finite element spaces of fc-forms. 

We require that the family of spaces X(T) fulfills the following consistency assumption: 

(4.2) tr r ,/ X(T) = tx T ',f X(T') whenever T,T' e T with / e A(T) H A(T'). 

In this case, we may define for any / G A(T), X(f) = try t X(T) where T G T is any simplex containing /. 
We also define X(f) as the subspace of X(f) consisting of all uj G X(f) such that tr f,dfu) — 0. Note that 

(4.3) tr g j X(g) - X{f) for all /, g G A(T) with / C g. 

Consequently, for each such / and g we may choose an extension operator Et g : X(f) — > X(g), i.e., a right 
inverse of tr g j : X(g) — > X(f). We say that a family of extension operators Ef g , defined for all /, g G A(T) 
with / C g, is consistent if 

(4.4) tr M = £/n 9 , s tr/,/n 9 for all /, g, h G A(T) with /, g C h. 
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In other words, we require that the diagram 



x(fng) 



X(h) 

tr 

X{g) 



commutes. 



One immediate implication of (4.4 1 is that for ui £ X(f), 
(4.5) tih, g Ef. h cu = Ef^ g uj for all /, g, h £ A(T) with / C g C h. 

A second implication is: 

Lemma 4.1. Let h £ A(T), and f,g £ A(h) with f <£. g. Then tth, g 0J = for all ui £ E f:h X(f). 



Proof. Let u> — Ef hfj, with /i £ X(f). Since / ^ g, we have / n g C 9/, and therefore tryj ng // = 0. Then, 
by p~4|, tT h>g uj = tr, iig -E/^M = E fn g ,g tr/,/n 3 M = 0. □ 



We now define an extension operator Ef : X(f) — » -Y(T) for each / 6 A(T). Given \x £ X(f), we define 
Ef/i piecewise: 



(4.6) 



Ef, T H if / C T, 
0, otherwise. 



We claim that for each g 6 A(T), tr g Ejfi is single- valued, so E/^i does indeed belong to X(T). To see this, 
we consider separately the cases / C g and / ^ g. In the former case, if T £ T is any simplex containing g, 
then / C T, and so 

tr T, s [(^/M)|T] = tr T , g E fiT n = E ft gn 



by (4.51. Thus tr<r, s [(i?//i)|T] does not depend on the choice of T containing g, so in this case we have 
established that ti g EfH is single- valued. On the other hand, if / < £ g then txT, g [{E '/ (i)\t] — f° r any T 
containing g, either because / ^ T and so (Ejh)\t — 0, or by Lemma 4.1 if / C T. Thus we have established 
that all traces of Effi are single- valued, and so we have defined extension operators Ef : X(f) — > X(T) for 
each / e A(7~). We refer to as the global extension operator determined by the consistent family of 
extension operators. 



We easily obtain this variant of Lemma 4. 1 



Lemma 4.2. Let f,g £ A(T), f <£ g. Then tr B u = /or a// w € E f X(f). 



Proof. Pick T £ T containing g. If / £ A(T), then we can apply Lemma 4.1 with h — T. Otherwise, 
(EfriW = for all fj, £ X(f). □ 



The following theorem is the main result of this section 
Theorem 4.3 



Let T be a simplicial triangulation and suppose that for each T 
subspace X(T) o/A fc (T) is given fulfilling the consistency assumption (4.2 I 

family of extensions operators E g j for all f,g£ A(T) with f C g. Define Ef, f £ A(T) by (4.6) 



£ T , a finite- dimensional 
Assume that there is a consistent 
Then the 



space X(T) defined in (4.1 1 admits the direct sum decomposition 



(4.7) 



X(T)= EfX(f). 
/eA(T) 
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Proof. To show that the sum is direct, we assume that X)/eA(T) u f = 0; where G Ef^X(f), and prove 
by induction that ajy = for all / G A(T) with dim/ < j. This is certainly true for j < k, (since then 
A fc (/) and, a fortiori, X(f) vanishes), so we assume it is true and must show that ui g — for g G A J+ i(T). 
By Lemma |4~2"1 

= tr f( H w /)= tr s w s- 
/eA(T) 

Hence, w g = £J g tr s w g = 0. We thus conclude that the sum is direct, and X(T) D (B/eA(T) E fX(f). 
To show that this is an equality, we write any u> G X(T) in the form 

n-1 
j=k 

where ui k — u>, and for k < j < n, u 3 G X(T) is defined recursively by 

/6A 3 (T) 

We shall prove by induction that for k < j < n 

(4.8) tr f uPeX(f), f G Aj(T). 

Assuming this momentarily, we get that oj 3+1 — lo 3 G X2/eA (T) E fX(f). Also, w"|t = tr*r<^™ € X(T) for 
all TeT, and w n = EtgT ti'T^lr)- Thus, ui G /eA (T) E f,rX(f) as desired. 



To prove (4.8 1 inductively, we first note it is certainly true if j = k, since A(/) = A(/) for / G Aj,(T). 
Now assume (4.8) and let g G A J+ i(T). We show that tY g ui 3+1 G X(g), by showing that tr^w :,+1 = for 



h G Aj(g). In fact, 



tr/, w^ 1 = tr/, w J - tr/, tr/ up . 



4.2 



tr/, Ef tif lo 3 = unless 



Now tTfUj 3 G X(f) by the inductive hypothesis, and therefore, by Lemma 
/ = h, in which case tr/, Ef tr/ lu 3 — tr/, a; 3 . Thus, 

tr/, oj 3 ^ 1 = tr/, — tr/, w J = 0. 

This completes the proof of the theorem. □ 



Remark. By considering the case of a mesh consisting of a single simplex T, we see that 

(4.9) X(T)= %!(/). 

/e a(t) 



The decomposition (4.7 1 is very important in practice. It leads immediately to a local basis for the large 
space X(T) consisting of elements Ef\i, where / ranges over A(T) and fi ranges over a basis for the space 

x(f). 

We close this section with the simplest example of this theory. Let X(T) = V r (T) — V r A°(T) be the 



polynomial space discussed in Section 2.3 Then (4.2 1 is fulfilled and the trace spaces X(f) are simply V r (f) 
for / G A(T). For f,g G A(T) with / C g, the trace operator tr ff j and barycentric extension operator Et g 
are given in barycentric coordinates as follows. If a G Nj] : ™ with [a] C 1(g), then 



tr s ,/(AT 



(X f ) a if 1(f), 
0, otherwise. 
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For a e Ng :n with \a\ = r and [a] C 1(f), then E ftg (X f ) a = (X 9 ) a . We now check that the family of 
barycentric extension operators is consistent, i.e., we verify ( |4.4| . We must show that if f,g,hd A(T) with 
f,gCh, then 

tr h , g E f , h (\f) a = E fngtg tr f , fng (\'r 
for all multi-indices a with \a\ — r and [a] C 1(f). Indeed, it is easy to check that both sides are equal (A 9 )" 
if [<*] C 1(g) and zero otherwise. Note that, in this case, the decomposition (4.9 1 is simply the Bernstein 



decomposition ( |2.6| . If we then define, as in the general definition (4.1) above, 

■p r A°(T) = { to e L 2 (fl) | uj\ T G V r (T) VT e T, tv f u is single-valued for / e A(T) }, 



then the decomposition (4.7) gives a decomposition of the space "P r A°(T), i.e., the space of continuous 
piecewise polynomials of degree < r. 



5. Degrees of freedom and the dual decomposition 



Although our main interest in this paper is obtaining direct sum decompositions for polynomial differ- 
ential forms that are analogous to the Bernstein decomposition for ordinary polynomials, we include here 
a discussion of another decomposition, referred to as the dual decomposition, for completeness and as an 
illustration of the general theory developed in the previous section. 

Before we consider the case of differential forms, we review the corresponding decomposition for poly- 
nomials. For the construction of finite element spaces based on the local space V r (T), a basis for the dual 
space V r (T)* is given, with each basis element associated to a subsimplex of T. This is referred to as a set 
of degrees of freedom for V r (T). The degrees of freedom then determine the interelement continuity imposed 
on the finite element space. Indeed, in the classical approach of Ciarlet [3J, the degrees of freedom and 
their association to subsimpliccs is used to define a finite clement space. For this purpose, what matters is 
not the particular basis of "P r (T)*, but rather the decomposition of this space into the spaces spanned by 
the basis elements associated to each simplex. For the standard Lagrange finite elements, this geometric 
decomposition of the dual space is 

(5.1) V r (T)* = W r (T,f), 

feA(T) 

where 

W r (g, f) := { V; e V r (g)* I <AM = J{tt g , f u)ri, 77 € 7Vdi m /-i(/) }• 
We note that for ui G V r (h), txhj w is uniquely determined by © g£ A(/) W r (h,g). 



Consequently, if / C h £ A(T), we may define an extension operator Ff t h — Fj h : V r (f) — > V r (h), 
determined by the conditions: 



(toh, g F fih uj)r) = j (tr /i9 cj)77, rj € V r - dimg -i(g), g E A(f), 
iP(F fjh Lj)=0, ^eW r (h,g), 9 eA(li), g£f. 

To apply the theory developed in Section[4] we need to check that the extension operator is consistent, i.e. 



that it satisfies (4.4 1. For f,g C h, let u G V r (f), and set ^ := Ffng,gtof,fngU & 'Pr(g), v := tr/j, 9 Ff^u € 
V r (g). For any face e C g D /, tr 9)e /i = tr/ i6 lu = tr g , e v. Therefore ip(fj-) = ^(v) for all "0 € W r (g, e) with 
e e A(</) such that e C /. Also, for e S A(<7) with e ^ /, it follows from the definition of the extension 
that for all tp € W r (g,e), ip(fi) = = ^(v). Thus we have shown that the extension operators Ff : h form a 
consistent family. The decomposition 

Vr(T)= F f , T [V r (f)], 

feA(T) 
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corresponding to (4.9 1, is now called the decomposition dual to (5.1l. Furthermore, from Theorem 4.3 we 



obtain a corresponding direct sum decomposition for the assembled space V r (T) = V r A°(T) of the form 



(4.7) 



In the remainder of this section, we present analogous results for the spaces V r A k (T) and V r A k (T). This 
will be based on the following decompositions of the dual spaces V r A k (T)* and V~ A k (T)* , established in 
[2 §4, Theorems 4.10 and 4.14]. 

Theorem 5.1. 1. For each f G A(T) define 

W k (T, /) := | TP G V r A k (Ty | = J tr TJ LuAvfor some r, G V; +k _ dim f A d ™ f- k (f) 

Then the obvious mapping 'P~ +k _ Ail - a fA dlm ? ~ k (f) — ► W k (T,f) is an isomorphism, and 

T r A k (T)*= W k (T,f). 
/e A(T) 

2. For each f G A(T) de/me 

W*~{T, /) := | V G V;A k {T)* | ^(w) = | tar,/ wAi) /or some v G P r+Jk _ dta/ _ 1 A dim '-*(/) 
TTien the obvious mapping Vr+k-dim /-iA dlm f~ k (/) — > W k ~ (T, f) is an isomorphism, and 



V-A k (T)*= W k ~(T, f). 

/6A(T) 

Note that as in the polynomial case, if uj G V r A k (T), then tr T , / cj is determined by the degrees of freedom 
in W k (T, g) for 5 G A(/). In particular, if uj G T J r A k {T) such that all the degrees of freedom associated to the 
subsimplices of T with dimension less than or equal to n — 1 vanish, then w G T J r A k {T). The corresponding 
property holds for the spaces P~A k (T) as well. 

An immediate consequence of this theorem are the following isomorphisms, that will be used in the 
following section. 



Corollary 5.2. 



V r A k (T)* = V; +k _ n A n - k (T) and P,7A fc (T)* = V r+k _ n _ 1 A n -"(T). 



As in the case of 0-forms, if / C h G A(T), we define an extension operator : T r A k (f) — ► V r A k (h), 
determined by the conditions: 

(tah, g Fj£u)r) = / (tr/, g w)??, 77 G Pr-dimg-lCg), 3 G A(/), 

9 Jg 

^»=0, i/>eW?{h,g), geA(h), g <£ f. 

We may similarly define an extension operator Ff'^ : T , ~A k (f) — > T'~A k (h). The verification of the 
consistency of these families of extension operators is essentially the same as for the space V r (T) given 
above, and so we do not repeat the proof. 
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6. BARYCENTRIC SPANNING SETS 



Let T — [xq, . . . ,x n ] C K™ be a nondegenerate n-simplex. The Bernstein basis described in Section 2.3 
above is given in terms of the barycentric coordinates {Ai}™ =0 C V\(T). The main purpose of this paper is 
to construct the generalization of the Bernstein basis for the polynomial spaces V r A k (T) and V~A k (T). In 
the present section, we will give spanning sets and bases for these spaces and for the corresponding spaces 
with vanishing trace expressed in barycentric coordinates. Note that the bases given in this section depend 
on the ordering of the vertices. These are not the bases we suggest for computation. 

For convenience we summarize the results of the section in the following theorem, referring not only to 
the n-dimensional simplex T, but, more generally, to any subsimplex / of T. Here, we use the notation 

(6.1) dX{ = d\i {1) A • • ■ A d\{ [k) G Alt fe T f 



for / G A(T), cr G E(l : k, : n) with [cr] C 1(f), and (j)£ for the Whitney form defined in ( [63] ). 
Theorem 6.1. Let f G A(T). 

1. Spanning set and basis for V r A k (f). The set 

{ (X f ) a d\l | a G \a\ = r, a G E(l : k, : n), [a, a] C 1(f) } 
is a spanning set for V r A k (f), and 

{ (X f ) a dXl | a G K n , \a\ = r, a G E(l : k, : n), [a, a] C 1(f), min[a] > minl(/) } 

is a basis. 

2. Spanning set and basis for P r A k (f). The set 

{ (Xf) a dXl | a G K n , M = r, a G E(l : k, : n), [a, a] = } 
is a spanning set for V r A k (f), and 

{(Xf) a dXl\aeN o Q : n ,\a\ = r,aeZ(l:k,0:n), la, cr] = 1(f) a, = if i < mm[l(f) \ {a}} } 
is a basis. 

3. Spanning set and basis for V~ A k (f) . The set 

{ (XfTti | a G K n , l«l = r - 1, a G E(0 : k, : n), [a, a] C Z(/) } 
is a spanning set for V~A k (f), and 

(6.2) {(AO Q ^|aGN° : ",|a|=r-l,crGS(0:fc,0:7i), [a, a] C 1(f) , a, = i/i<min[cr]} 

is a 6asis. 

4- Spanning set and basis for V~ A k (f) . The set 

{ (X f ) a 4>l | a G K n , M = r - 1, a G E(0 : fc, : n), [a, a] = Z(/) } 
is a spanning set for V~ A k (f) , and 

{ (X f ) a 4 | a G |a| = r — I, a G E(0 : fc,0 : n), [a, a] = Z(/), a.^Oi/K min[a] } 
is a basis. 
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6.1. Barycentric spanning set and basis for V r A k (T). Observe that dXi G Alt 1 E". Furthermore, 
d\i(xj — y) = Sij for any y in the subsimplex opposite Xj. In particular, tr^j dXi = for any subsimplex 
/ G A(T) with Xj ^ / or equivalently i G I(f*). Furthermore, {dXi}™ =0 is a spanning set for Alt 1 ^™, and 
any subset of n elements is a basis. Therefore, writing dX a for dX^, the set 

{dX a | cr G E(l : k,0 : n)} 

is a spanning set for Alt fc M" , and the set 

{ dX a | cr G E(l : fc, 1 : n) } 
is a basis. The forms dX a G Alt fc R™ have the property that for any / G A(T), 

tr Tj/ d\(j = if and only if [cr] n J(/*) ^ 0. 
More generally, for polynomial forms of the form \ a d\ a G V r A (T), we observe that 

tv TJ {X a d\ a ) = if and only if [a,(r]nl(f) ^ 0. 
In particular, X a dX a G V r A k (T) if and only if [[a, cr] = {0, . . . , n}. 



Taking the tensor product of the Bernstein basis for V r (T), given by (2.3), with the spanning set and 
basis given above for Alt W 1 , we get that 

Proposition 6.2. The set 



{ \ a d\ a | a G N° - n , |a| = r, a G E(l : k, : n) } 
is a spanning set for V r A k (T), and 

{ X a d\ a | a G N° Q :n , \a\ = r, a G E(l : jfc, 1 : n) } 

is a 6asis. 

Restricting to a face / G A(T), we obtain the spanning set and basis for r P r A k (f) given in the first part 
of Theorem 16.11 



6.2. Barycentric spanning set and basis for V r A k (T). For / G A(T) and a G E(0 : k, : n) with 
[cr] C define the associated Whitney form by 

fc 

(6.3) ^ = EMr^W < (0) A • • • A dAi w A • • • A dA^ (fc) . 

i=0 

Just as we usually write Xi rather than A^" when the simplex is clear from context, we will usually write <j> a 
instead of (p^. We note that if k = 0, so that the associated subsimplex f a consists of a single point Xi, then 
4> a = Xi. It is evident that the Whitney forms belong to V\A k (T). In fact, they belong to A k (T). This 
is a direct consequence of the identity 

(6.4) KdXp — (0) 5 



which can be easily established by induction on k, using the Leibniz rule (3.4 1. In fact, the set 

{ <t> a | a G E(0 : k, : n) } 

is a basis for A k (T) . Furthermore, tTxj <t>a — dX a ^ A • • • A dX^^ is a nonvanishing constant fc-form on 
/ = f a , while tr T j ($> a - - for / G Aj-(T), f =/= f a - Therefore, we refer to </> CT as the Whitney form associated 
to the face f„ . 

For cr G E(0 : k, : n) and < j < k, we let <p a j be the Whitney form corresponding to the subsimplex of 
f a obtained by removing the vertex o~(j). Hence, 

j-i 

0<rj = y^-l^A^rfA^p) A • • • A dX a (i) A • • • A dX a( j) A . . . dA CT(fe) 
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- ^2 (-^-W) Ar(0) A • • • A d\ aU) A • • • A d\ a{i) A... d\ a[k) . 

From this expression, we easily obtain the identity 

k 

(6.5) Y^rtfKwfas = g G E(0 : fc, : n). 

i=o 

Correspondingly, for j ^ \o\ , we define 

^j'ct = Xjd\ a — dAj A <f> a . 

Thus, modulo a possible factor of —1, 4>j a is the Whitney form associated to the simplex [xj, f a ]. For these 



functions, we obtain 



(6.6) 



J2 fa* = ( x ^ dX ° - ( Z A ^ 



= ( x j) d ^ + ( dA j) A ^ = (Z A J-) dA ^ = dX °- 

Now consider functions of the form A a <fi a , where a G Nq", \a\ = r — 1, <j G E(0 : fc, : n). It follows from 
the relation (3.5) that these functions belong to V~A k {T). In fact, they span. From the identity (6.5) we 



know that these forms are not, in general, linearly independent. The following lemma, cf. [2, Lemma 4.2], 
enables us to extract a basis. 

Lemma 6.3. Let x be a vertex of T. Then the Whitney forms corresponding to the k-subsimplices that 
contain x are linearly independent over the ring of polynomials V(T). 

Using these results, we are able to prove: 
Proposition 6.4. The set 

(6.7) { X a c/) a | a G N° Q - n , \a\ = r - l,a G E(0 : A,Q : n) }. 
is a spanning set for V~ A k (T) , and 

(6.8) { X a 4>a- 1 a G N° - n , \a\=r-l, a G E(0 : fc,0 : n), a, = £/i < wxx\a\ } 
is a basis. 



Proof. Let a G N$ n , with |a| = r - 1, and p G E(0 : fc - 1, : n). The identity (6.6) implies that 



\ a d\ 



,. — Z X '°J!'- 



and hence all forms in V r -xA k (T) are in the span of the set given by (6.7). Furthermore, if a G E(0 : fc, : n), 



we obtain from (6.4 1 that 



K (\ a dX a ) + \ a <p a (0) = A Q </v, 

and therefore all of n[P r -iA k+1 (T)] is also in the span. By the definition of the space r P~A k (T), it follows 
that (6.7 1 is a spanning set. To show that (6.8 1 is a basis, we use the identity (6.5 1 to see that any form 



given in the span of (6.7 1 is in the span of the forms in (6.8 1. Then we use Lemma 6.3 combined with a 



simple inductive argument, to show that the elements of the asserted basis are linearly independent. For 
details, see the proof of Theorem 4.4 of [2]. □ 



Restricting to a face / G A(T), we obtain the spanning set and basis for V r A k (f) given in the third part 
of Theorem 16.11 
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6.3. Spaces of vanishing trace. In this subsection, we will derive spanning sets and bases for the corre- 
sponding spaces of zero trace. This will be based on the results obtained above and Corollary |5.2| which 
leads to the dimension of these spaces. We first characterize the space P~A k (T). 

Proposition 6.5. The set 

{ \ a (f> a | a e N° :I \ \a\ = r - 1, a € S(0 : fc, : n), [a, aj = {0, . . . , n} } 
is a spanning set for V~ K k (T) and 

{ A Q 0(j | a G Nq : ™, \a\=r-l, a G E(0 : fc, : n), [a, aj = {0, . . . , n}, a, = if i < min[ff] } 
is a basis. 



Proof. Since [a, cr] = {0, . . . , n), each of the forms X a c/> rT is contained in V r A k (T). Moreover, the condition 



aj = if i < minjer] reduces to er(0) = in this case. Lemma 6.3 implies that the forms X a (f> <T for which 



<r(0) = are linearly independent. The cardinality of this set is equal to (?) dimV r - n +k-i which is equal 
to dimV~ A k (T) by Corollary 5.2 This completes the proof. □ 



Restricting to a face / € A(T), we obtain the spanning set and basis for V r A fc (/) given in the fourth 



part of Theorem 6.1 



Finally, we obtain a characterization of the space V r h k {T). 
Proposition 6.6. The set 

{ X a d\ a | a e N° :n , |a| = r, a € E(l : k, : n), [a, aj = {0, . . . , n} } 
is a spanning set for V r K k (T), and 

(6.9) { \ a d\ a | a G N[! : ", \a\ = r, a G E(l : fc,0 : n), [a, a] = {0,...,n}, a* = i/i < minlcr*] } 
is a &asis. 



Proof. Since [oj,erfl_= {0, . . . , n}, each of the forms X a dX a is contained in V r A k (T). Furthermore, we have 
seen in Corollary g that V r K k {T) S 7V +fc _»A"~ fc (r), whence, dimP r A fc (T) = (^)( r + fc ). On the other 
hand, the cardinality of the set given by (6.9 1 can be computed as Y] ^ Aj ■ Bj, where Aj is the number of 
elements a G E(l : fc,0 : n) with minjcr*] = j, and for each fixed such cr, Bj is the number of multi-indices 
a satisfying the conditions of (6.9), namely 

Wrfi and B,-( r+ ;:j- r 



fc - j 

Hence, the cardinality of the set is given by 

k 

E 



3=0 



n-j 
k-j 



r + k — j — 1 
n-j 



r - 1 
n — fc 

r- 1 

71 — fc 

Here the first identity follows from a binomial identity of the form 

:t)Q- (:)(::: 



E 

J'=0 

r+k 



r + fc — j — 1 
r - 1 



while the second is a standard summation formula. Hence, the cardinality of the set given by (6.9 1 is equal 



to the dimension of P r A k (T). To complete the proof, we show that the elements of the set (6.9 1 are linearly 
independent. Denote the index set by 

S := { (a, a) G N[j :n x E(l : k, : n) \ \a\ = r, {a, a} = {0, . . . , n}, a, = if i < min[(r*] }, 
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so we must show that if 

(6.10) c aa \ a d\ a = 0, 

for some real coefficients c QCT , then all the coefficients vanish. Since the Bernstein monomials A" are linearly 



independent, (6.10 1 implies that for each a £ Nq with \a 



(6.11) 



E 



dX„ = 0. 



First consider a multi-index a with a$ > 0. Then the definition of the index set S implies that minjer*] = 
for all the summands in (6.111. Since the corresponding dX a are linearly independent, we conclude that all 
the c aiJ vanish when ao > 0. Next consider a with ctQ = but a\ > 0. If (a, a) £ S, then minjer*] = 1, 
and again we conclude that c aa = 0. Continuing in this way we find that all the c aa vanish, completing the 
proof. □ 



Restricting to a face / £ A(T), we obtain the spanning set and basis for V r A k (f) given in the second 
part of Theorem |6.1| 



7. A GEOMETRIC DECOMPOSITION OF V~A k (T) 

In this section, we will apply the theory developed in Section [4] with X(T) = V~A k (T) to obtain a 
geometric decomposition of r P r A k (T) into subspaces Ef[V~ A k (f)}, where Ej is the global extension operator 
constructed as in Section [4] from a consistent family of easily computable extension operators. The resulting 
decomposition reduces to the Bernstein decomposition (2.6 1 in the case k = 0. 



We first note that if T,T' £T with / £ A(T) n A(T') then tr T>/ V~A k (T) = tr T ',/ p-A k (T') = V~A k {f). 

k,r.— 



Hence, the assumption (4.2) holds. Furthermore, for f,g G A(T) with / C we define E — E^ g 

r-A k (f) 

(7.1) 



V r A (g) as the barycentric extension: 



(\f) a <pl » (xm, Officii/). 



This generalizes to k- forms, the barycentric extension operator Ej T on V r , introduced in Section 2.3 Since 



the forms (A^)"0£ are not linearly independent, it is not clear that (7.1l well-defines E. We show this in 
the following theorem. 



Theorem 7.1. Th 



ere is a unique mapping 



E = Ef£~ from V-A k (f) to V-A k (g) satisfying (fl \ 



Proof. We first recall from part 3 of Theorem |6.1| that the set 

{ {\ f ) a (f>l | a £ N° :n , \a\ = r - 1, a £ E(0 : Jfe,0 : n), [a, a] C 1(f), a, = if i < mm{aj } 



is a basis for V~A k (f). Hence, we can uniquely define an extension E by (7.1 1, if we restrict to these basis 
functions. We now show that (7.1 ) holds for all (\f) a (j>l with [a, a] C 1(f). To see this, we use the identity 
(6.5). Consider forms (\^) a <p } a which do not belong to the given basis, i.e., s := minja] < minjo - ]. Write 
A" = A^ and let p £ E(0 : k + 1,0 : n) be determined by [p] = {s} U [cr]. Then, by dgj), 



fc+i 



(\r^ = Y^i-iy^x^x^s, 



and so 



fc+i 



fc+i 



(x f rti = ^(-ly-'ixffx^^, (xtk = E(- 1 ) J '" 1 ( A9 )^ A "(i)^- 
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Hence 

fc+i fe+i 

3=1 3=1 

and the proof is completed. □ 

Theorem 7.2. The family of extension operators E is consistent, i.e., for all f,g,h G A(T) with f,gQh, 
and all uj € V~A k (f), 

tr/i, 9 Ef^w — Ef n g. g tifjng UJ - 

Proof. It is enough to establish this result for uj — (\f) a (j>l, with [a, cr] C 1(f), since such cj span r P~A k (f). 
Now for such pairs (a, a), E f>h [(\ f ) a <l>l} = (X h ) a (j>^ and then 



tt h , g E frh [(\f) 



if la,aj Cl(/nj), 
0, otherwise. 



On the other hand, 



and hence 



tr 



/./riff 



(X f ) a 4>l = i( xfn9 ) a ^ fn9 )<r, if [a, ff ] Cl(/n S ), 

1 0, otherwise, 



p tr rA A^/ J( A9 ) Q ^' if[a >0 -]cj(/n^), 

*Wr/,,n s (A') Q otherw . s ^ 



□ 



From Theorem 4.3 we obtain the desired geometric decomposition of 7^ A fe (T). 
Theorem 7.3. 

P r -A fc (T)= E f [V-A\f)]. 



/eA(T) 

dim />fc 



where Ef : V r A k (f) — > V T A k (T) denotes the global extension operator determined by the family E k ' g 



The final part of Theorem 6.1 furnishes an explicit spanning set and basis for V~A k (f), and so this 
theorem gives an explicit spanning set and basis for V^A k {T). We discuss these explicit representations 
further in Section [9l 

We now turn to a geometric characterization of the extension operator E : T J ~A k (f) — > V~A k {T, /). 
To this end, we say that a smooth fc-form u> g A k (T) vanishes to order r at a point x if the function 
x <— > uj x (v\, . . . , Vk) vanishes to order r at x for all v\, . . . , Vk € K", and that it vanishes to order r on a set 
g if it vanishes to order r at each point of the set. Note that the extension operator E = E*jT has the 
property that for any /i E r P 7 TA k (f), E k 'J^~ fi vanishes to order r on /*. In fact, if we set 

V~A k (T, f) = {ioe V~A k (T) | uj vanishes to order r on /* }, 

we can prove 

Theorem 7.4. V~A k (T,f) = E[P~A k (f)] and for fi e V~A k (f), Efi = E)'Jj~ fi can be characterized as 
the unique extension of fi to V^A [T, /). 
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Proof. We note that the second statement of the theorem follows from the first, since tixj from E[P~ A k (/)] 
to P~A k (f) has a unique right inverse. Since E[P~A k (f)] C P~A k (T, /), we only need to prove the opposite 
inclusion. Without loss of generality we may assume that / = [x m+ i, . . . ,x n ], f* = [xq, . . . ,x m ], for some 
< m < n. We proceed by induction on m. When m = 0, we may assume without loss of generality that 
the vertex x is a * the origin. Now V~ A k = V r -\A + K?i r _ 1 A fc+1 (where H. r denotes the homogeneous 
polynomials of the degree r). Since KH. r -iA k+1 C Ti r A k , every element lo € nH. r -iA k+1 vanishes to order 
r at the origin. On the other hand, no non-zero clement of V r -\A k vanishes to order r at the origin. Thus 
V-A k (T,f) = KH r -iA k+1 . It follows from Theorem 3.3] that 

dimV-A k (T,f) = dim K n r ^A k+1 (T) = fclj) ( ? + \ ~ ^ = ^K^(f) = dimE[p-A k (f)], 

and since E[P~A k (f)] C p-A k {T,f), E[P~A k {f)] = p-A k (TJ). 

Now suppose that lo vanishes to order r on the m-dimensional face [xq, ... ,x m \ with m > 0. Let T' — 
[xi, . . . ,x n ), lo' — ttTTtLO. Then lo' 6 V~A (T') vanishes to order r on the (m — l)-dimensional face 
[xi, . . . ,x m ], so, by induction, lo' = Ef^'t^ for some /i G P,rA fc (/). Furthermore, since w vanishes to order 
r at xo, we can use the result established above for m = to conclude that lo — Et'.t^' — Et> .tE f ^r 1 li- 



However, it follows immediately from (7.1 I that Ej» j*Ef y = Et j., and hence the two spaces are equal, and 



the theorem is established. □ 



8. A GEOMETRIC DECOMPOSITION OF V r A k (T) 

In this section, we again apply the theory developed in Section [ij this time with X(T) — V r A k (T). In 



this case condition (4.2) is obvious, since tr^ / V r A k (T) = V r A k {f). In view of the previous section, one 
might hope that we could define the extension operator as 

{X f ) a dX{ ^ X a dX a . 

However, this does not lead to a well-defined operator. To appreciate the problem, consider the space 
7- , 2A 1 (T), where T C R 2 is a triangle spanned by the vertices xq, xi,X2, and let / = [xi, x-i\. Then \[\2(d\{ + 
dX f 2 ) = 0, but AiA 2 (dAi + dX 2 ) = -X 1 X 2 dX Q ^ 0. 

To remedy this situation, we will show that for f,g € A(T) with /Cj,a consistent extension operator 
E = Ef g : V r A k {f) - V r A k {g) is given by 

(8.1) (X f ) a d\l ~ (ATC /,S , I a > CT 1 ^ 

where <tp£>f> 9 is defined as follows. We first introduce forms ip" ' 9 £ Alt 1 T g defined by 



(8.2) ^J,9 =dx i-^ d\°, iei(f), 



and then define V"' 7 ' 9 G Alt fe T g by 

(8.3) W J ' B =€(i)*''-W%k)> ^€E(l:fe,0:n), [a]Cl(/), 

A geometric interpretation of ip^>f<9 will be given below. 

First we show that i? is well-defined and is, in fact, an extension operator. 
Theorem 8.1. There is a unique mapping E — E/ r g : V r A k (f) — > V r A k (g) satisfying (8.1). Moreover it is 



an extension operator: tl g j Ej'^lo — lo for lo G V r A k (f) 
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Proof. By the first part of Theorem |6.1[ the set 

{(X f ) a dX f a \a£ N° - n ,\a\ = r, a £ E(l : k, : n), [a, a] CI(/), min[(r] > minX(/) } 
is a basis for P r A fe (/) . Hence, we can define an extension E by (8.1), if we restrict to the basis functions. We 



now show that (8.1l holds for all (X^) a dXl with [a, <r] C 1(f), i.e., also when minjfr] = minX(/). Writing 
dX a = dA CT (i) A dX p , and using the fact that X^ez(/) = on the face /, we can write 

dX f a = - dX f j AdX f p . 

Hence, 

c,f,g a ,/,a,/,S _ f\g\a-La,f,g 



E{(Xf) a dXl} = -E{(Xf) a ]T dX f J AdXj,] = -(XT £ M*? 1 '" = WW, 

jex(f) jei(f) 

where in the last step we have used the fact that X^ex(/) ip?'^ 9 = 0. 
That 23 is an extension operator follows directly from the observation 

which holds since X^ez(/) dXj = on the face /. □ 

Theorem 8.2. The family of extension operators E is consistent, i.e., for all f,g,h £ A(T) with f,gQh 
and all uj £ V r A k (f), 

Proof. It is enough to establish this result for u> — (A^) q g?A£, with [a, u] C X(/). Now for such pairs 
(a, a), E s . h [{X f ) a {dX f ) a ] = (X h ) a ^ J ' h . To determine tx h ,g[WWv f ' % we consider three cases. When 
[a] C J(>), but [a] £ X( 5 ), tr^pT^'M = 0, since tr M [(A' l ) Q ] - 0. If [a] C J(/ n <?), then 
t r /i,g[(A' l ) Q ] = (A 9 )", so we need only compute trh, g ^"'^' h - We do this by first considering tr^ jS ip"'^' h for 
i £ 1(f). If i £ l(f n gi), we have 

tr hi9 ^' h =tr h , fl (dA?-^ £ dA£) = dA?-0 £ dA* = 

On the other hand, if i £ 1(f) \X(fr\ g), then since a$ = 0, ^p"'^' h = dX^ and so tr/j, g t/j"'^' h = 0. Combining 
these results, we obtain 

10, otherwise. 

But 

V ns )°dA£ nfl > if [a,*} CX(fng), 



tv fJng [(Xf) a dXi] = | 
and hence 

^/n 9 , ff tr /Jng [(A0 Q ^] 



0, otherwise, 

= |(A 9 ) Q C /n9 ' 9 , if[a,o-]CJ(/n<7), 
] 0, otherwise. 



□ 



From Theorem 4.3 we obtain the desired geometric decomposition P r A (T). 
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Theorem 8.3. 

V r A k (T) = E f [P r A k (f)]. 



feA(T) 

dim />fc 



where Ef : P r h k (f) — > P r K k (T) denotes the global extension operator determined by the family E^ 



k,r 

g' 

Combining this result with the second part of Theorem |6.1[ we obtain an explicit spanning set and basis 
for 7VA fe (T) (see Section^. 

We now turn to a geometric characterization of the extension operator E = E/ g : P r K k (f) — > P r A fc (T, /). 
First, we will motivate the choice of E, and in particular the forms ip"^' 9 , by establishing some some 
additional properties of these forms. Observe that any multi-index a determines a convex combination of 
the vertices Xi of T, namely 

x a . — |o;| ^ ^ Oi m x m G T, 

m 

and if [a] C 1(f), then x a G /. For each such multi-index a, we then define the vectors 

t a i = x a - x t = S~] a m (x m -xi), l£l(f*). 
\a\ 

Clearly, for each such a, W decomposes as the direct sum Tf © span{f () ; 1 1 6 where denotes the 

tangent space of /. See Figure 8.1 This decomposition defines a projection operator P = Pf tCt : R™ — * re- 



determined by the equations Pv — v for v G Tf and Pt a i =0 for I G X(f*). Hence, we have 
(8.4) P* f a A\t k T f = {ae Alt fe R n \ aj t at = 0, 2 G 1(f)}. 

Furthermore, since d\j(x m — x{) — Sj m for any j G 1(f), rn G 1(f), and I G 1(f*), we get for [a] C 1(f), 



.5) i>i ,i,T (totl) = d\ l (t ai ) - pr E (*<*') = ~~ ( c 



iexCf) 



It follows that 

(8.6) ^ f ' T (v ir -- ,v k )=^ f > T (P Vl ,--- ,Pv k )=ti T j^' f ' T (Pvi,--- ,Pv k )=d\ f a (Pvi,--- ,Pv k ). 
Hence, in the language of pullbacks, 

W f ' T = P}, a d\l, 

where Pf a is the pullback of P/ jQ , , and so 

E Lg [(X^rdXi] = (XTP*f, a d4- 




Figure 8.1. T = [x ,x 1 ,x 2 ,x 3 ], /= [x ,xx], a= (3,1,0,0), R 3 = T f © span{i 2ct , t 3a }. 
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Recall that the geometric characterization in the previous section hinged upon the fact that a form in 
V~A k (T) which vanishes to order r on /* and has vanishing trace on / must vanish identically. Now this 
is not true for an arbitrary element of the larger space V r A k (T). Returning to the example given at the 
beginning of this section, where T = [xq, x\, x%] and / = [351,2:2], the form 

lo = X 1 X 2 [dX 1 + dX 2 ] = -X 1 X 2 dX e V 2 A l (T) 

vanishes to second order at /* = {xq}. However, trxjco also vanishes. Thus, additional conditions on to will 
be needed in order to insure that lo is uniquely determined by tiT,f to. We say that to vanishes to order r + 
on /*, if a; vanishes to order r on /* and the following conditions hold: 

(8.7) d%oojt al = 0, I G 1(f), H C 1(f), \a\ = r. 

Here := Iljei(/) ^t-i with dt jl = tji ■ V, the directional derivative along the vector tji :— Xj — X\. The 
contraction operator j is defined at the start of Section |2.4| 

Note that d tjl \ = hj for i,j G 1(f), j G l(f*). It follows that if a,(3 G NJ, with |a| = \0\ = r, 
[a],[/3]CI(/), J and?eX(.D,then 

(8.8) 9^A Q = fora^/3and9°A a = a!. 
Setting 

V r A k (T, f) = { uj G V r A k (T) | lo vanishes to order r + on /* }, 
we can now give the geometric description of the extension operator E. 

Theorem 8.4. V r A k (T,f) = E[V r A k (f)] and for \x G V r A k (f), Efi = Ef^fj, can be characterized as the 
unique extension of [i to V r A k (T, /). 



Proof. We note that the second statement of the theorem follows from the first, since tirj from E[P r A k (f)] to 
P r A k (f) has a unique right inverse. To prove the first statement, we first show that E[P r A k (f)] C V r A k (T, f). 
Observe first that E[(X^) c "(dX^) a ] — X a ip%'f' T vanishes to order r on /* since (A) Q does. Next, note that 
8 ) tells us that c*f [X a ip"'f' T ] = a!-0"'^' T if /? = a and vanishes if f3 is any other multi-index of order r with 



] Q 1(f)- Therefore, the conditions (8.7 1 for vanishing of order r + are reduced to verifying the conditions 
ip£'f' T _,t a i = for all / G l(f*). However, this follows immediately from the definition of the wedge product 



and (8.51 



To show that V r A k (TJ) C E[P r A k (f)}, we use Lemma 
admits a representation of the form 

00= ^ 

Hci(/) 

\a\—r 

for some a a G Alt k I 



2.1 



to see that any element lo G V r A k (T 1 f) 



81 and 



However, invoking 
on /*, then a a jt a i = for all I 6 T(f*), and hence by 
6]), that co G E[V r A k (f)]. 



7), we conclude that, if lo vanishes to the order r + 

It therefore follows from 
□ 



4|, a a G P; Q 



Alt fe T 



9. Construction of bases 



From Theorem |7.3| (4.6 1, Theorem |7.1[ and part 4 of Theorem 6.1 one immediately obtains explicit 
formulas for a spanning set and basis for V~ A k (T), with each spanning and basis form associated to a 
particular face / G A(T). The forms associated to / vanish on simpliccs T 6 T that do not contain /, while 
for T containing /, the spanning and basis forms are given by 

{ (X T ) a ^ | a G K n , \a\ = r - 1, a G S(0 : k, : n), [«, a] = J(/) } 
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and 

{ {\ T ) a 4>l I a e N° : ", \a\ = r - l,a € S(0 : k, : n), [a, a] = !(/), a ( = if i < min[<r] }, 

respectively. Note that the spanning set is independent of the ordering of the vertices, while our choice of 
basis depends on the ordering of the vertices. Other choices of basis are possible as well, but there is no one 
canonical choice. 



The same considerations give an explicit spanning set and basis for V r A k (T), based on Theorem 



8.3 



(|4.6|), Theorem 8.1 and part 2 of Theorem 6.1 The corresponding formulas for the spanning set and basis 

{ (A r rC /,T I a £ K n , M = r, a € E(l : k,0 : n), 



are 



and 



{ (A T rC /,T I « G K n , H = r, a e E(l : k, : n), [a, <r] = J(/), a, = if i < min[/(/) \ [cr] ] }, 



respectively, where ip"'^ T is defined by (8.2) and (8.3 1. 

Bases for the spaces V~A k and V r A k are summarized in Tables 9.1 - 9.4 for n = 2, 3, < k < n, r = 1, 2, 3. 
In the tables, we assume i < j < k < I, and recall that the Whitney forms and <^jjfc are given by: 



— Aj g^Aj — Aj d\i 



4>ijk = Ai d\j A ciAfc — Aj d\i A dA& + A/c d\i A dAj 



Table 9.1. Bases for the spaces TV A 1 and T^A 1 , 





TV A 1 






r 


Edge [xi,Xj] Triangle [xi,Xj,x k ] 


Edge [a^ay] 


Triangle [xi, xj , x k ] 


1 


4>ij 


A j c£A j , Aj d\-i 




2 


X k <f>ij, ~Xj4>ik 


X 2 dX,, A 2 dA, 
Xi\jd(Xj — Ai) 


XiXjdX kl XiXfcdXj 
XjX^dXi 


3 


{A 2 , A 2 , Xi^j}<j>ij {Ai, \j,\k}\k<t>ij 
{Ai, Xj , X k }Xj<j> ik 


AfdV,-, A^A, 
XfX j d{2X :j - X t ) 
XiXjdiXj - 2Ai) 


{Ai, Aj ■, AfclAiAjdAfc, 
{Ai, Aj ;, Xk}XiXkdXj 
{Aj , Afc}Aj X^dXi 



Table 9.2. Bases for the spaces V~ A 1 and V~A 2 , n = 3. 



TV A 1 



TV A 2 



r Edge [xi,Xj] Face [xi, xj ,x k ] Tet [x^ Xj, x*,, x;] 



Face [xi, Xj, xj,] 



Tet [xi, Xj , Xfc, xj] 



{Ai, Aj}0i 



{A 2 ,A 2 ,A,Aj}^j 



Afe</>ij, Aj</>i fc 

{Ai, Aj, Afc}Afc0ij 
{Ai, Aj, AfcjAj-c/iij; 



AfcA;0ij 

XjXi<j>ik 

^j^k4>il 



<t>ijk 

{Ai, Aj, X k }4>ij k 



Ai^ijfc, X k 4>iji 

^j&ikl 



{A 2 , A 2 , A|}</)ij fc {Ai, Aj, A fe , Ai}A ; 0ij fc 

{A,;Aj ;, AiAfc, Aj Xk}cj>ijk {Ai, Aj, X k , Xi}X k rj>iji 
{Ai, Aj , X k , X t }Xj(j>iki 
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r Edge [x{, Xj] Face [xi, Xj , x k ] Tet [xt, Xj , x k , xi] 

1 XidXj , \jd\i 

2 X 2 dXj, \ 2 -d\i, Xi\jd(Xj — Xi) XiXjdX k , XiX k dXj, XjX k dXi 

3 X^dXj, X?Xjd(2Xj — Xi) {Xi, Xj}XiXjdX k , XiXjX k d(2X k — Xi — Xj) XiXjX k dXi, XiXjXidX k 
X^dXi, XiX?d(Xj — 2Xi) {Xi, Xk}XiXkdXj , XiXjX k d(2Xj — Xi — X k ) XiX k XidXj, XjX k X[dXi 

{Xj, X k }XjX k dXi 



Table 9.4. Basis for the space V r A 2 , n = 3. 



Face [xi,Xj,x k ] Tet [xi,Xj , x k , x t ] 



1 X k dXi A dXj , XjdXiAdXfc, XidXj A dX k 

2 X^dXi A dXj , XjXfrdXi A d(X k — Xj) X k X[dXiAdXj, XjX[dXiAdX k 
X 2 dXiAdX k , XiXjd(Xj — Xi) A dX k XjX k dXi A dXi, XiXidXjAdX k 
X^dXj A dX k , XiX^dXj A d(X k — Aj) XiX k dXj A dXi, XiXjdX k A dXi 

3 X 3 k dXiAdXj, X 3 -dXiAdX k , XjdXj AdX k {X k , A;}A fc A ; cZA; A dXj 
X 2 X k dX t A d(2X k - Xj), XjX^dXi A d(X k — 2Xj) {Xj, X k , Xi}XjXidXi A dX k 
XfXjd(2Xj - Xt) AdX k , X 2 X k dXj A d(2X k - A,) {Xj, X k , Xi}XjX k dXi A dXi 
XiX 2 d{X 3 - 2Xi) A dX k , XiX 2 dXj A d(X k - 2A;) {Xi, Xj, X k , AJAjA^A., A dX k 

XiXjX k d(2Xj — Xi — X k ) A d(2X k — Xi — Xj) {Xi, Xj, X k , Xi}XiX k dXj A dXi 



{Xi, Xj , X k , Xi}XiXjdX k A dXi 
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